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Finding the Complete Path and Weight Enumerators 

of Convolutional Codes 

I. Onyszchuk 1 

Communications Systems Research Section 

A method for obtaining the complete path enumerator T(D, L, I ) of a convo- 
lutional code is described. A system of algebraic equations is solved, using a new 
algorithm for computing determinants, to obtain T(D, L, I) for the (7,1/2) NASA • 

standard code. Generating functions, derived from T(D,L,I), are used to upper 
bound Viterbi decoder error rates. This technique is currently feasible for constraint j 

length I\ <^10 codes. . . j 

A practical / fast algorithm is presented for computing the leading nonzero 
coefficients of the generating functions used to bound the performance of constraint 
length I< < 20 codes. Code profiles with about 50 nonzero coefficients are obtained 
with this algorithm for the experimental I< = 15, rate 1/4 , code in the Galileo 
mission and for the proposed K — 15, rate 1/6, *2~dB code. 

As an alternative, algorithms have also been developed to 
calculate the first few coefficients (the distance and bit er- 
ror profiles) of these enigmatic generating functions [4,5]. 
Unfortunately, lists of these numbers are bulky (see Ta- 
bles 2 and 3), and the minimum number of terms required 
to approximate the decoder error bounds depends upon 
the code rate and channel noise level. 

A code’s complete path enumerator T(D,LJ) con- 
tains the number of paths having identical triples: weight, 
length, number of input Is. The least-magnitude pole of 
the weight enumerator T(D) = T(D, 1, 1) determines the 
point at which the union bounds [1] diverge, while addi- 
tional poles and residues yield the dominant terms in the 


I. Introduction 

Convolutional codes such as the (7,1/2) NASA stan- 
dard have been used for satellite and deep-space commu- 
nications during the past 20 years. In 1971, Viterbi [1] de- 
fined generating functions for upper bounding error proba- 
bilities of convolutional codes on memoryless channels. In 
practice, for codes with more than eight states, these func- 
tions are still unknown, so error bounds have been eval- 
uated using numerical matrix multiplications [2,8], which 
require extensive computations for each channel noise level. 


1 Also a student in the Electrical Engineering Department, Califor- 
nia Institute of Technology, Pasadena, California. 
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By Cramer’s Rule, 


partial fraction expansion of T(D). Concise, yet very accu- 
rate, analytic approximations to several generating func- 
tions’ coefficients may be obtained from these dominant 
terms. Furthermore, the poles and residues of T(D) may 
eventually help unlock the structure of convolutional codes. 

The generating functions T(D) and T(D, L, 7) for the 
(7,1/2) NASA standard code (among many others) were 
obtained with a simple determinant algorithm described in 
this article. In addition, a new algorithm is presented for 
computing the initial coefficients of T(D) y dT(D , L, I)/dL 
and dT(D, L, I)/dI at L — I— L These numbers are used 
to upper bound a Viterbi decoder’s event, node, bit, and 
symbol error rates. The techniques explained here may 
be adapted, with modifications for code nonlinearity and 
Euclidean instead of Hamming distances, to find weight 
enumerators of trellis codes [3,5]. 


II. The Complete Path Enumerator 

A binary, rate k/n , feed-forward convolutional en- 
coder is defined by kn binary generator polynomials 
gij(x) y each representing the transfer function from the 
*th input to the j th output (x is a delay operator). Let 
g_ be the vector whose rth component is the coeffi- 
cient of x r in gij(x). The encoder’s memory is m — 
ZlLi [deg 9ij( x )]i which for rate 1/n codes is K — 1. 
If the zth encoder memory cell contains S{ then the en- 
coder is in state s = YT= The encoder’s state 

diagram is a directed graph whose edges (corresponding 
to branches in the associated trellis diagram) have labels 
n,- j = D d LI h if there is an edge from state j into state i 
while a,j = 0 otherwise. During a transition from state 
j into state z, b is the number of ones input to the en- 
coder and the number of ones output by the encoder is 
d — ]T£ =1 g_ lh • [2 j + (f mod 2)] 2 , where * is a modulo-2 
inner product and [k] 2 denotes the binary representation 
of the integer k. 

Let X $ be the trivariate generating function of all sim- 
ple paths: those from state 0 into state s via nonzero states. 
Xq counts all simple paths into state 0, called fundamental 
paths. Note that X s is indexed by the destination state 
(s) while the source state (0) is constant. Now define A 
as the 2 m — 1 by 2 m ~ 1 adjacency matrix of the encoder 
graph with state 0 and its edges removed. The entry in 
row i and column j of A is a,-j, which is nonzero only if 
there is a directed edge connecting state j to state i (they 
are adjacent). Now the following set of linear equations is 
constructed: 

A [X\ , X 2i X 2 m - i] T = [&i, 0 i a> 2 ,o* a 2 m -i,o] 


x _ det(Aj) 

* det(A) 

for i > 0, where A ; is A with all column i entries a r> * 
replaced by a r>Q for all rows r. A code’s complete path 
enumerator ( by weight, length, and number of input Is) 

2 fc 

T(D,L,I) = X 0 = J2a 0ttj X, j 
3 = 1 


where {zj} are the 2 l states having edges into state 0 (so 

flO.Zj ^ 0 ). 

As an example with k = 1 for simplicity, the m = 2 
encoder in Fig. 1 is in state s =2si + s 0 . The corre- 
sponding state diagram in Fig. 2 leads to the equations 

A[X 1 ,X 2) X 3 ] T = [a li o,0,0] T : 


1 —LI 0 ' 


’Xi 


~ D 2 Ll~ 

-DL 1 -DL 


x 2 


0 

1 

C3 

1 

o 

C5 

1 

1 


/ 3j 


0 


Now T(D, L, I) = X 0 - D 2 LX 2 = D 5 L 3 I/(l — DLI — 
DL 2 I), so the code’s weight enumerator is 

n5 00 00 

r(D)= TT2D = Y. 2d ~ SDd = E P«O d 

d=5 d=d,„. 

where p(d) = 2 d “ 5 is the number of weight d fundamental 
paths, and d fre€ = 5 is the code’s free distance. 


Ill- Reducing an Adjacency Matrix 

When k = 1, there are only 3(2 m -l)*-2 nonzero out 
of (2™ — l) 2 entries in A,» and they are located in a special 
pattern. Since A, is sparse if m > 3, reduction is useful be- 
fore computing the determinant. The following example, 
different from the one in the previous section, illustrates 
the reduction procedure. The code with generator polyno- 
mials gn(D) = 1 + D + D 3 and gi 2 (D) = 1 + D + D 2 + D 3 
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has an adjacency matrix A (with L = I = 1 to simplify 
entries) 


(B 2 ".-i) and B as the resulting upper left 2 m 2 by 2 m 2 
submatrices, 


1 

0 

0 

(D 2 )- 1 

0 

0 

0 


1 

0 

0 

-1 

0 

0 

-1 

0 

1 

0 

-D 2 

0 

0 

0 

-D 

0 

(0)1 

0 

-D 

0 

0 

-D 

0 

0 

1 

-D 

0 

0 

0 

-D 

0 

0 

1 

—D 

0 

0 

-D 

0 

0 

0 

1 -D 


The round brackets (parentheses) above indicate values in 
Aa^-i that are different from those in A. Since k = 1, 
the determinants of only these two matrices are needed for 
T(D,L,I ), and this notation will lead to their simultane- 
ous computation. 

As in Gaussian elimination, -a r [ r / 2 j times row |r/2j 
is added to rows r = 2 to 2 m -l so that A becomes zeroed 
below its main diagonal for columns 1 to 2 m_1 — 1 : 



i 

0 

0 

m -i 

0 

0 

0 



0 

1 

0 

(D 4 ) -D 2 

-1 

0 

0 



0 

0 

1 

(D 2 ) -1 

—D 2 

0 

0 



0 

0 

0 

(D 5 )l-D 3 

-D 

-D 

0 



0 

0 

0 

(D 5 ) -D 3 

1 -D 

-D 

0 



0 

0 

0 

(D 3 ) -D 

-D 3 

1 

—D 



0 

0 

0 

(D 3 ) -D 

-D 3 

0 

1 -D 


Therefore, 

det (A) equals the determinant of 

the resulting 


lower right 2 m_1 by 2 m_1 submatrix. To further reduce A, 
row |r/2j times — a r ,2’ n - 1 +L r / 2 J ' s ^ded to each row r = 
2 m — 3 to 2 m_1 (4 to 5 here), so that columns 2 m_1 +2 m_2 
to 2 m — 2 (both 6 here) are zeroed above the main diagonal. 
Now det (A) = det (A), where A 2 m-i) and A are the new 
lower right 2 m_1 by 2 m_1 submatrices: 


(£> 4 + D 5 )l- 

-D 2 -D 3 

- D-D 4 0 

-D 2 ’ 

(D 4 + D 5 )~ 

-D 2 — D 3 

1 -D-D 4 0 

-D 2 

(£» 3 ) 

-D 

-D 3 1 

-D 

(D 3 ) 

-D 

-D 3 0 

1 -D 


To zero column 2 m_1 of A (above the diagonal entry 
a 2 , »- 1 , 2 m -' = a 2 m -i, 2 m -i 7^ 0)> — 5r i 2 r "- 1 /a2 m_I ,2 m_1 times 
row 2 m_1 is added to each row r = 1 to 2 m_1 — 1. Define 


(D 4 + D 5 -D 6 )l-D-D 2 -D 3 + D 4 -D+D 2 -D 4 
(D 4 + D 5 -D 6 ) — D 2 — D 3 +D 4 1-2 D+D 2 -D 4 

1 - D 


This reduction method simultaneously produces two dense 
2 m-2 by 2 m_2 matrices (I3 2 „,-i ), B. with the same determi- 
nants as the corresponding original sparse 2 m 1 by 2 — 1 

matrices (A^™-*)) .A. 


IV. A Determinant Algorithm 

The following algorithm yields the determinant of any 
N x N matrix B having entries from a Euclidean domain, 
such as the set of all polynomials with integer coefficients. 
A sequence of matrices B 1 N _ 1 \ B' N 2 \ . . ., B 11 1 is 

computed with each matrix having the same determinant 
(up to sign) as B. Starting with B (N) = B, step j in the 
algorithm produces the numerator of each entry in 

row i and column k of B^ -1 ^: 


for j = N to 2 


(step) 


for i = j— 1 to 1 (row index) 

for k — j to 1 (column index) 


b U) b U) _ b U ) 6 (i) 


Naturally, if 6 j+ + i!i-M = 0 P rior to ste P then an y 
column k such that ^ 0 must first be interchanged 

with column j- hi. This operation negates the determinant, 
so a counter t is incremented to record the event. If no 
such column k exists, row j + 1 is zero, so the algorithm 
is stopped and det (i?) = det ( B ^) — 0. Also note that 

^Ar-M,Ar+l = ^ initiall y* 

The following example with N = 4 illustrates the 
above reduction procedure. When j ’ = 4, after 

MV * an<1 M? times row 4 are added to rows 

3, 2, and 1, respectively (corresponding to j~ N = A in the 
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algorithm), in order to zero all entries in column 4 above 
kW 

°44 > 


B (3 ) = 


»8> 

A(3) 

°12 

6 ( 3) 

°13 


& 

*S? 

A(3) 

°2X 



6 (4) 

°44 

b( 4 ) 
0 44 

b (4) 

°44 

A( 3 ) 
°31 

°32 

b (3) 

°33 

b {4) 

°44 

b {4) 

°44 

b (4) 

°44 

b (4) 

U 41 

b {4) 

°42 

b^ 

°43 


& (4) 

°44 


Then, after the j = 3 step during which —b^/b^ and 
~b{ 3 ^ /b^ times row 3 are added to rows 2 and 1, 


B (2) = 


€ 

6 ( 2 ) 

°2l 

b (3) 

°33 

fc(3) 

°31 

6 ( 4 ) 

°44 


a(<) 


6£ 

« 

a(2) 

°22 

6 (3) 

°33 

6 (3) 

p 32 

6 f4) 

°44 


"41 


42 




& (4) 

°44 


Since only elementary row operations on B have been 
performed, det(B) = det(B w ) =s det(B (2 ^). The final 
step ( j — 2) produces a B ^ matrix identical to b ( 2) ex- 
cept with b^/b^ 0 0 0 in the first row. Therefore, 

det(B) = det(B^) = b[\\ This result generalizes to any 
nonsingular N x N matrix B . 

Lemma. Entry b[ l \ = det(B^) = ( — 1)* det(B), 
where t is the total number of column interchanges per- 
formed by the algorithm in computing 6^. 

Proof. For each value of j from N to 2, when k ~j, 
the algorithm makes = 0 for all 1 < i < j - 1 so 

that all column j entries above the main diagonal become 
0. Since B^ is zero above its main diagonal, det(B) is 
the product of (— 1)* and the diagonal entries in b/ 1 ^. 

det(B) = (-1)' det(B (1 >) 


= (- 1 )' b[]l • 

Expanding b\ 3 J b ( f] - b^J using the equation in 

the algorithm shows that divides this expression so 

there are never any remainders. This is expected because 
the algorithm implements a recursive factorization of the 
determinant written as a sum of products of matrix entries. 
The algorithm differs from standard Gaussian elimination 
because the diagonal entries in the reduced matrices are 
not made equal to 1. Also, after calculations with a par- 
ticular value of j are completed, all entries’ denominators 
are previous pivot numerators. 


V. Path Distance, Length, and Bit Error 
Approximations 

The complete path enumerator T(D, L, 1) for the 
m = 6, rate 1/2, NASA standard code was obtained by 
using the preceding algorithm to simultaneously compute 
determinants of the two 16 by 16 reduced matrices 
and B_. The 76 poles of T(Z)) for this code are plotted on 
the complex plane (Fig. 3) along with the unit circle for 
reference. Using only the six least-magnitude poles (indi- 
cated by the large points in Fig. 3), an approximation to 
the partial fraction expansion of T(D) is 


T{D) 




+ 


-f 


^2 


ocdD 1 -j- ajD 1 — 


+ : 


T2 


+ 


l + a b D 1 -a* b D 1 + aJDj 

00 

= D 10 [ r i ( a d)* + ri(-a d ) k + r 2 (o t )* 


*=o 




+ r a (-a»)* + r2(a;)* + r;(-a;) 


D 


OO 

= [2ri(a d ) 2i +4Re{r 2 (a J ) 2 *}] £> 2i+I ° 
i= 0 

where a d = 2.3876225 is the reciprocal of the least- 
magnitude pole of T(D) t oc b = 1.657193e“° 983418 > is the re- 
ciprocal of the pole with next smallest magnitude, * means 
complex conjugate, and j — yj — 1 . The residues are 


N 


= (-o' n 


s«> 

jj 


J-J- M+ 1) 

i=i °j+ij+i 


- p i a l) 

n = d ’ 


CIa 




and 


p _ -^K 1 ) 

2 ^Q'K 1 ) 
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where D l0 P(D)/Q(D) = T(D) and Q\D) is the deriva- 
tive of Q(D). The other poles in Fig. 3 may be ignored 
because their magnitudes are greater than 0.8, and the 
corresponding residues have magnitudes less than 0.07. 

Define p(d),i(d) y and £(d) as the coefficients of D d in 
T{D),dT(D,L,l)/dL, and dT{D } LJ)/dL , respectively, 
at L = I = 1 (these generating functions are shown in 
the Appendix). The number of fundamental paths hav- 
ing weight 2k + 10 is 

p(2*+10) « 6.82(2.3876225) 2 * 


For the additive white Gaussian noise (AWGN) channel 
with bit signal-to-noise ratio E^/Nq, 

P d = Q(y/2dRE b /N 0 ) 

where Q(x) is the Gaussian integral function [2, p. 62]. On 
a binary symmetric channel with crossover probability p, 


(2i — 1\ x2*-l-e ^ /2i — 1 

p 2 , = p^-i = £ ( e )p e (i-p) <( i 

e=i ' ' 


+ 4.25(1. 657193) 2t cos (0.310- 1.967/:) 

Similarly, the terms corresponding to the six least- 
magnitude poles of T(D) in the partial fraction expan- 
sions of dT(D, L y I)/dL and <9T(D, L, I)/dI at L = J=1 
were used to obtain the approximations 

£(2k+ 10) « (77.725 +22.625fc)(2.3876225) 2 * 

+ 39.3(1. 657193) 2 *cos (0.485- 1.967/:) 

+ (2*+l)7.2676(1.657193) 2i cos (0.383 - 1.967*) 

t(2* + 10) « (24.474 12.018Ar)(2.3876225) 2 * 

+ 9.942(1.657193) 2 *cos (0.575- 1.967*) 

+ (2*+ 1)2. 8723(1. 657193) 2i cos (0.366 - 1.967*) 

which have a relative error < 0.0001 for k > 4. 

A rate 1/n Viterbi decoder’s bit error rate (BER) on 
a binary-input, output-symmetric [2], discrete memoryless 
channel is bounded by 

oo 

BER < £ i{d)P d 

d = di rC€ 

where P d is the probability that the decoder outputs a fun- 
damental path having distance d from the one transmitted. 
The probability that a 6-bit symbol is decoded incorrectly 
is bounded by [9] 

SER& < £ K 6 - 1 - m M<0 + *(<01 p d 

d = d ftct 


[6]. For decoders using integer metrics, as for example on 
a binary-input, output-quantized AWGN channel, P $ can 
also be computed exactly [2, p. 291]. 

VI. Algorithms for Profiles of Convolutional 
Codes 

Finding the complete path enumerator of codes with 
memory greater than 8 currently seems infeasible. In these 
cases, an algorithm for distance profiles [4,5] may be used 
to calculate the first few nonzero coefficients of the gen- 
erating functions used for error bounds. However, these 
methods are fairly complex and some require extra compu- 
tation to ensure that the output is correct. In this section, 
Viterbi’s algorithm, with survivors replaced by vectors of 
integers that count paths, lengths, or bits, is applied on 
a noiseless channel to compute p(d) ) £(d) i and b(d) values. 
Rate 1/n codes are treated first to simplify the discus- 
sion. Define out 0 [s] and outi[s] as the number of ones 
that the encoder outputs going from state sO = [s/2 J and 
si = sO -h 2 m_1 into state s. Analogous to a state met- 
ric, the entry in row s > 0 and column 0 of a matrix W , 
referred to as W[$][0], will contain the least weight of any 
simple path with length < T trellis branches. For t = 1 
to coeffs (a parameter described later), W[s][<] will be the 
number of simple paths of weight W[s][0]+< —1 and length 
< T branches into state s > 0. For state 0, W[0][<] is 
always kept at 0, except W[0][1] = 1. The entries in a 
second matrix B count either the total number of ones 
input to the encoder (when the variable len = 0) or the to- 
tal length in trellis branches (when len = 1), of all simple 
paths having length < T (again B[0][t] = 0 always). These 
matrices are obtained for successive values of T starting 
with 1 by extending, one branch length at a time (an al- 
gorithm ‘step’), the code trellis starting from state 0 only. 
Thus the longest pathlength (T) explored by the algorithm 
equals the number of ‘steps’ executed. The algorithm ter- 
minates after step T * when W has reached values that 
will never change, which also forces B to remain constant. 
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Then since W[s][0] is the least weight of any simple path 
into state s > 0, d free = W[2 m ~ l ][0] 4 outi[0]. Also, 

W [2 m - l )[d--dte e +l] = p(d) 

5[2 m - 1 ][rf-i free + 1] = b(d) (£(d) if len = 1) 
for d — t/frcc to ~f* coeffs — 1 , 

Two versions of the basic algorithm above are pre- 
sented in C language format in Algorithm 1 and Table 1. 
In Algorithm 1, matrices P and A store previous W and 
B entries corresponding to simple paths of length < T, 
which are used to compute new W and B matrices for 
length < T-j-l simple paths. When change remains 0_af- 
ter step T* f W (and thus B) will never change because 
W[s][<] = P[s][i] for all s and t. .P[s][0] is initialized to 999 
for $ > 0, .P[0][1] = 1, and all other array values are ini- 
tialized to 0. If any second array index t + offset is < 0 in 
the W[s][2] and £[$][£] instructions, the array referenced is 
simply ignored. 

The algorithm requires storage for 2 m+1 (coeffs -j- 2) 
integers and the amount of work per step is proportional 
to this number. The number of steps executed, T * , equals 
the length, in trellis branches, of the longest fundamental 
path(s) having weight d free + coeffs - 1. The parameter 
coeffs should be set equal to [lOn/fc] because using this 
many nonzero terms in the union bounds gives results with 
three significant digits of precision when the bounds are 
tight enough to be useful. 

Setting coeffs = 0 and ignoring offsetO, offset 1, A, 
and B yields a simple and fast algorithm for finding df ree . 
About 2 m+1 bytes of storage and 2.1 CPU seconds (on 
a computer executing 12 million instructions per second) 
were required to obtain d free for the m = 14 Galileo and 
“2-dB” codes. 

For rate k/n codes with k > 1, a state s is parti- 
tioned into si, S2 j * • • j $k where s* corresponds to the con- 
tents of the ith encoder shift register. The output weight 
as the encoder enters state s is outj[s] and the input is 


j G [0 ... 2* — 1]. New Wr[«][<] and fl [«][*] values are com- 
puted using at most 2 k entries from each of the P and A 
matrices. 

For codes with 2 k « 2 m (such as rate 1/n with 
m > 3), Algorithm 1 may be improved by looping through 
groups of 2 k states called butterflies (see Fig. 4) instead 
of individual states and by computing W[s][^] and £[s][t] 
in place [7]. This reduces the storage memory required by 
almost one-half because the “double-buffering” matrices P 
and A are eliminated. If the entire vector W[2 m ’ _1 ] remains 
the same for m consecutive steps, it will never change be- 
cause there is a trellis path of length < m branches between 
any two states. The algorithm shown in Table 1, which 
incorporates these improvements, produced the first lOn 
nonzero coeffs in the profiles of two K — 15 codes: the 
rate 1/4 Galileo code (Table 2) and the proposed “2-dB”, 
rate 1/6 code (Table 3). These profiles took only a few 
minutes of CPU to generate and required storage of 7 and 
8.5 Mbytes, respectively, when 4 bytes were used for each 
integer. These memory requirements could be further re- 
duced by storing each integer in the smallest number of 
bytes needed (1 for W[s][0] to W[s][15] in the K = 15 codes 
above). 


Algorithm 1. A simple profile algorithm 


do { 

for (s = 1 to 2 m — 1) { 

50 = [5/2J ; 5l = 50 + 2 m_I ; 

bit = len + (1 —len) * (5 mod 2) ; 

W[s}[0] = min (P[s0][0] + outofa], P[5l][0] -f- outofs]); 
offsetO = W[s][Q] - P[s0][0] - out 0 [s] ; 
offsetl = W[a][0] - P[sl][0] - outo [*]; 
for (f = 1 to coeffs) { 

W[s][t] = P[*0][t + offsetO] -j- P[5l][t + offsetl]; 
B[s][i] = v4[^0][t -f offsetO] + 

>i[5l][t -f offsetl] -f bit * VF[s][f]; } 

} change = 0 ; 
for (5 = 1 to 2 m — 1) 
for (£ = 0 to coeffs) 

if (P[s][t] # VFfsfft]) { change = 1; 

PMM = WHM; B[s][t] = A[s}[t}; } 

} while (change ^ 0); 
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Table 1. A profile algorithm for rate 1 fn codes 


k — 0; stop = 0; 
do { 

for (£ = 0 to coeffs) 

{ temp[t] = W[2 rn ~ l ~ k [t]\ £[0][<] = W[0][t] = 0; } 

W[0][l] = l; 

for (s = 0 to 2 m_1 - 1) { 

sO = 5 » k\ (cyclically) 
si = sO -f 2 m_1 ~ k 

t*m = min (W[s0][0] + out 0 [2e], W[sl][0] + out! [2s]); 

twl[0] = min (iy[s0][0] -f outi[2« + l], W[sl][0] + outo[2* + l]); 

offsetOO = fu/0[0] - lV[s0][0] - outo[2s]; 

offsetOl = fu/0[Q] — W[«l][0] — out\[2s\\ 

offset 10 = fu/l[0] - W[jO][0] - outi [2s-j-l]; 

offset 11 = tud[0] - VT[sl][0] - outo [2a+'ljj 

for (< = 1 to coeffs) { 

tu;0[t] = V^[s0][t + offsetOO] + W[si][t + offsetlO]; 

twl[t] = W[sQ][t + offsetlO] + W[sl][t + offset 11]; 

tb0[t] - B[sO][t + offsetOO] + B[sl][t + offsetlO] + len * tw0[t]-, 

f6l[t] = B[s0][< -f offsetlO] + B[sl][t + offsetll] + twl[t]; 

for (* = 0 to coeffs) { 

W[sO][t] = twO[t]; W[sl][t] = twl[t]; 

B[sQ][t] = tb0[t]; B[*l][t] = <M[t] ; } 
k = k+ 1; if (k = m) k = 0; change = 1; 
if (iy[2 m -l-fc][ coe ffs] > 0) change = 0; 
for (t = 0 to coeffs) 

if (iy[2 m “ 1 “*][<] ^ temp[t]) change = 1; 
if (change =1) stop = 0; 

if (change = 0 and stop < m-1) {change = 1; stop+ + ; } 

} while (change ^ 0); 



Table 3. Rate 1/6 "2-dB" code profiles 


Table 2. Galileo code profiles 


Distance 

d 

Fundamental 
paths p(d) 

Bit errors 
»■(<*) 

Total 

lengths £(d) 

35 

2 

6 

7 

36 

1 

2 

5 

37 

4 

16 

22 

38 

2 

8 

12 

39 

3 

11 

17 

40 

5 

20 

28 

41 

6 

24 

46 

42 

17 

76 

122 

43 

24 

126 

214 

44 

29 

180 

285 

45 

39 

255 

438 

46 

66 

416 

721 

47 

94 

628 

1071 

48 

121 

850 

1478 

49 

175 

1313 

2260 

50 

277 

2086 

3643 

51 

415 

3361 

5855 

52 

639 

5304 

9388 

53 

934 

8010 

14161 

54 

1273 

11452 

20271 

55 

1906 

17550 

31381 

56 

2878 

27332 

49172 

57 

4054 

39750 

71705 

58 

5978 

60788 

109808 

59 

8864 

92738 

167861 

60 

12966 

139556 

253134 

61 

18984 

210112 

383008 

62 

27949 

317798 

581467 

63 

41092 

479512 

878975 

64 

60126 

720858 

1323152 

65 

87799 

1080933 

1987235 

66 

128712 

1622990 

2992979 

67 

189880 

2451782 

4530508 

68 

278589 

3682496 

6817868 

69 

408780 

5534126 

10261968 

70 

598271 

8283100 

15386816 

71 

875283 

12380669 

23050515 

72 

1286052 

18596544 

34662286 

73 

1888299 

27885609 

52045238 

74 

2768375 

41727376 

78013493 

75 

4057688 

62421220 

116865844 

76 

5953416 

93419654 

175122289 

77 

8732134 

139709066 

262220198 

78 

12809968 

208928290 

392628663 

79 

18786484 

312181796 

587384902 

80 

27548175 

466271448 

878292728 

81 

40412499 

696477455 

1313354906 

82 

59269748 

1039725314 

1962719710 


Distance 

d 

Fundamental 
paths p(d) 

Bit errors 

i(d) 

Total 

lengths 1(d) 

56,57 

1,5 

2,15 

3,19 

58,60 

1,3 

2,12 

2,14 

61,62 

5,12 

25,56 

35,84 

63,64 

11,5 

43,24 

67,40 

65,66 

8,11 

44,62 

68,95 

67,68 

8,11 

48,62 

76,98 

69 

27 

167 

267 

70 

30 

162 

277 

71 

36 

216 

363 

72 

54 

366 

573 

73 

74 

464 

785 

74 

89 

610 

998 

75 

94 

670 

1104 

76 

126 

912 

1524 

77 

163 

1209 

2022 

78 

226 

1676 

2814 

79 

290 

2236 

3785 

80 

369 

2920 

4993 

81 

493 

4051 

6846 

82 

574 

4780 

8168 

83 

767 

6571 

11236 

84 

979 

8562 

14687 

85 

1182 

10474 

18250 

86 

1574 

14282 

24860 

87 

1996 

18516 

32193 

88 

2618 

24594 

43183 

89 

3407 

32955 

57577 

90 

4238. 

41914 

73499 

91 

5353 

53757 

94399 

92 

7006 

71430 

126401 

93 

8932 

92712 

164631 

94 

11418 

120946 

214330 

95 

14401 

155175 

275986 

96 

18467 

202902 

361135 

97 

24039 

268439 

479664 

98 

30325 

344146 

616671 

99 

38662 

446878 

800288 

100 

49690 

583672 

1048171 

101 

63930 

762130 

1371587 

102 

81742 

990268 

1785532 

103 

103839 

1278325 

2308219 

104 

133335 

1666564 

3012971 

105 

170357 

2159215 

3912282 

106 

217467 

2801764 

5081111 

107 

278512 

3640320 

6613934 

108 

356223 

4721974 

8592622 

109 

456347 

6135943 

11180051 

110 

583546 

7956498 

14517787 

111 

746528 

10327464 

18861798 

112 

954389 

13376948 

24468075 

113 

1220261 

17333391 

31748900 

114 

1562164 

22493842 

41242505 

115 

1997088 

29126250 

53477032 
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Appendix 

Generating Functions for the (7,1/2) NASA Code 

Using the determinant algorithm in section IV, the complete path enumerator T(D,L,I) for the m = 6, rate 1/2 
NASA code was found to contain 1529 numerator and 2799 denominator trivariate terms. The code’s weight enumerator, 

T(D) = T(D, 1, 1), is 

11D 10 — 6Z? 12 -25D 14 +£> 16 +93D 18 -157? 20 -176D 22 - 76 D 24 +243 D 26 +417Z? 28 - 228Z) 30 -1156 D 32 -49L> 34 
+2795 D 3 * +61 ID 38 -5841D 40 -1094D 42 +9575D 44 +1097Z? 46 -119007? 48 -678 D 50 +11218 D 52 +235L> 54 -80687? 56 
— 18£> 58 +4429D 60 -20D 62 -1838D 64 +8D 66 +562£> 68 -D 70 -120 D 72 +16D 76 - D 80 


1 _ 4 £)2 _e D 4 -30 D 6 +40D 8 +85D 10 -81D 12 -345D 14 +262D 16 +844D 18 -403£> 20 -1601D 22 +267D 24 +2 5 09 D 28 
+389 Z? 28 -3064 D 30 -2751 D 32 +2807D 34 +8344D 36 -1960 D 38 -16133D 40 +1184D 42 +21746D 44 -782D 46 -21403D 4 
+561I> 50 +15763D 52 — 331£> 54 - 87 66Z? 56 +131D 58 +3662D 60 -30Z? 62 -1123D 64 +3£> 66 +240£> 68 - 32 D 72 +2D 76 

- ll£)i° + 38D 12 + 193£> 14 + 1331D 16 + 7275D 18 + 40406Z? 20 + • • • 

The other two generating functions used to compute error bounds, dT(D, L, I)/dL and dT(D, L, I)/dI at L — I— 1, 
both have denominators equal to the square of T(D)’s denominator above. Their numerators are, respectively, 

121D 10 -387D 12 -706Z? 14 +1460Z? 16 +3970Z? 18 -6157D 20 -11643Z? 22 +8725D 24 +28677D 26 +12195 Z? 28 
+88D 30 -170654D 32 -306124D 34 +817895D 36 + 1637616Z? 38 -2879440Z? 40 -6106837Z? 42 +8568521Z? 44 +18636083Z) 4 
-22431469Z? 48 -48921504D 50 +52678351Z? 52 +113105887Z? 54 -112260733D 56 -232580537D 58 + 2 1 7 3 3 7 1 70D 
+426400859 D 62 -379787502Z? 64 -696667758D 66 +592954735D 68 +1013294336D 70 - 815739185D 72 -131112 4 721D 
+968225450D 76 + 1509511967D 78 -955561 827Z? 80 -1548537967D 82 +721812022D 84 +1419185285D 86 -302418615Z? 
-1166004400Z? 90 -173248817Z? 92 +861869027D 94 +545751792D 96 -574515412D 98 -713158178D 100 +345328990Z? ° 2 
+676371 119Z? 104 -186475599Z? 106 -515274530Z? 108 +89761092D 110 +326707300D 112 -38067910Z? 114 -174942675D 
+ 14010022D 118 +79516060D 120 -4391424D 122 -30654965D 124 H-1145504/} 126 -f 9969013D 128 -241287D 

-2706667D 132 +39344D 134 +603670Z? 136 -4651Z? 138 -107908Z? 140 +354Z? 142 +14883Z? 144 -13D 146 -1488D 
+96D 152 -3D 156 


36£>io -77 D 12 -140Z? 14 +813Z? 16 +269Z? 18 -4414Z? 20 +321Z? 22 +14884Z? 24 -5273D 26 -40509Z? 28 +39344Z? 39 
+83884 D 32 -177469D 34 -111029D 36 +608702 D 38 -29527D 40 -1820723Z? 42 +817086D 44 +4951082Z? 46 -3436675Z? 48 
— 12279246D 50 +10300306D 52 +27735007D 54 - 25648025D 56 -56773811D 58 +55659125Z) 60 + 104376199D 

— 106695512D 64 -170819460D 66 +180836818Z? 68 +247565043Z? 70 -270555690D 72 -317381295Z? 74 +356994415Z? 

+360595622 Z? 78 -41540 1723Z? 80 -364292177Z? 82 +426295756D 84 +328382391Z? 86 - 385686727D 88 -264812337D 

+307287819Z? 92 +191225378D 94 -215144035Z? 96 -123515898Z? 98 +131946573D 100 +71124860D 102 -70570661D 94 

-36310569D 106 +32722089D 108 +16308558D 110 -13052172D 112 -6380604 £> 114 +4433332D 116 +21475657? 
-1265046D 120 — 612040D 122 +297721D 124 +144665D 126 -56305D 128 -27569 D 130 +8232D 132 +40667? 134 - 874D 136 
— 435D 138 +60D 140 +30D 142 -2D 144 -D 146 
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